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Abstract: The precision of measurements of the strong couphng constant using event- 
shape variables in e"'"e~ annihilations is limited by theoretical systematic uncertainties. 
The uncertainties are related to missing higher orders in the perturbative predictions for 
the event-shape distributions. A new method is presented for the assessment of theoretical 
uncertainties in ag ■ This method evaluates the systematic uncertainty of the parameter Os 
from the uncertainty of the prediction for the distributions from which it is extracted. The 
perturbative uncertainties are calculated on a purely theoretical basis, without accessing 
measured distributions. The method is therefore especially suited for an unbiased com- 
bination of results from different observables or experiments. It is universal and can be 
applied to other processes like jet production in deep-inelastic ep scattering or in hadron 
collisions. 
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1. Introduction 

Studies of Quantum Chromodynamics in e'''e~ annihilations have been carried out over 
more than 30 years at increasing centre-of-mass energies. Event-shape variables have proven 
to be key observables in both annihilation and deep-inelastic scattering processes. The 
understanding of perturbative and non-perturbative aspects of QCD has grown with the 
study of event shapes. While event shapes were first introduced to characterise global 
properties, it was soon realised that their distributions are sensitive to the strong coupling 
constant ctg. The perturbative prediction for a generic infrared-collinear (IRC) safe event- 
shape variable y can be computed to second order in Og, 



with coefficient functions A and B obtained from integration of the ERT [1] matrix el- 
ements. Using this type of prediction, first determinations of as were performed at the 
PEP and PETRA colliders. In the era of the LEP experiments the calculations for cer- 
tain classes of variables were improved by resumming leading and next-to-leading loga- 
rithmic terms (NLL). These calculations, matched to fixed-order expressions, enlarged the 
kinematic range of applicability for ag extractions and reduced the systematic theoreti- 
cal uncertainty. Recently, event-shape variables have also been used extensively to study 
non-perturbative power law corrections. It has been shown that hadronisation corrections. 




(1.1) 
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scaling with inverse powers of the momentum transfer Q can be modelled with one or a 
few non-perturbative parameters [2]. These parameters can in turn be related to moments 
of an effective coupling at low scales, and have been extracted from the data [3]. 

At present measurements of Ug using event-shape variables and 0(a^)+NLL resummed 
calculations are available at centre-of-mass energies between 22 GeV and 206 GeV. For 
illustration a selection of these measurements is shown in Fig. 1. 
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Figure 1: Measurements of as from event-shape variables as function of centre-of-mass energy. 

New measurements re-analysing data from the JADE collaboration [4] were performed 
using a technique similar to that applied for a combined analysis by the LEP collaborations 
[5]. In particular, the same theoretical framework was used and at most of the energies 
a similar set of variables was combined. At almost all energies perturbative theoretical 
uncertainties dominate other systematic uncertainties, e.g. aX, Q = Mz the perturbative 
uncertainty is 5% and the quadratic sum of statistical, experimental and hadronisation 
uncertainties is 1.3%. 

It is the objective of this note to examine ways of estimating the theoretical uncer- 
tainties and to propose a comprehensive and consistent method to combine different un- 
certainty estimates. In addition, the proposed method should allow the calculation of the 
uncertainty associated with a given measurement without using the data from which ag 
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is extracted. This is notably important for global combinations of various measurements, 
where a coherent procedure must be applied. 

The origin of perturbative theoretical uncertainties is the truncation of the perturbative 
series at a given order of the coupling constant, whereby the truncated prediction acquires 
a dependence on the renormalisation scale. 

To illustrate this point, a fully inclusive observable, the total hadronic cross section 
(t(c^c~ qq), is considered as an example for which third order 0{ag) calculations are 
available. The ratio Rz of hadronic to leptonic widths at the Z boson resonance reads 
(with renormalisation scale ji = Q) [6] 



r(Z^ — > leptons) V TT / 

^ ' n=0 

with 

i?EW = 19.934 , Co = 1 , ci = 1.045 , C2 = 0.94 , C3 = -15 , (1.3) 

One estimate of the theoretical uncertainty on Rz is based on the size of the last term 
in the series, taking the difference between the full next-to-next-to-leading order (NNLO) 
prediction and the next-to-leading order (NLO) truncation. Taking the preliminary com- 
bined LEP measurement Rz = 20.767 ± 0.025 [7] the NNLO result is as{Mz) = 0.1240 ± 
0.0037(exp) and the central value shifts down to af^'^(Mz) = 0.1214 when the NLO 
prediction is used. Hence, the theoretical uncertainty estimated by the NNLO to NLO 
difference is Aa^ = ±0.0026. The 'true' error would be the difference between the NNLO 
result and the complete theory. The NNLO to NLO difference can only be a good estimate 
of that if the convergence of the perturbation scries is fast enough, i.e. if the size of the 
expansion coefficients remains of order unity. In the case of i?z, however, the size of the 
third order term amounts to 60 % of the second order term size. 

Of course the size of the last known term of the series is not the only way of esti- 
mating the uncertainty on the prediction. Another method originates from the fact that 
dimensional regularisation, used to define formally divergent loop corrections prior to MS 
renormalisation, requires the introduction of a renormalisation scale n at which the cou- 
pling, Qs(/i^), is defined. While the choice of scale is arbitrary, higher-order perturbative 
corrections contain terms proportional to powers of In^u^/Q^, which order-by-order com- 
pensate for the scale dependence of the coupling. Only in the presence of all orders this 
compensation is complete. For a calculation to order a" there is residual scale-dependence 
of order a""^^; since it would be cancelled by higher order terms, its size can be taken as 
indicative of the magnitude of missing higher-order contributions. 

The scale dependence of Eq. (1.2) is given by the dependence of the NLO and NNLO 
coefficients on the scale parameter = ^ as follows 

C2{x^) = cl + (i^ciXxixl (1.4) 
cz{x^) = 4- (2/3ocO + ci/3i) hixl + cil3l In^ , (1.5) 
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Figure 2: The dependence on the renormalisation scale parameter of as{Mz) using Rz at NNLO 
(solid line) and at NLO (dashed line). 



where the coefficients c° are those of Eq. (1.3) and the first three coefficients of the P 
function are 
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(1.6) 



The number of active flavours Uf is taken to be five at Q = Mz- The renormalisation 
scale dependence of the coupling constant can be parameterised in the modified minimal 
subtraction scheme at 3-loop level [8] as a function of the scale A = A^^^ 



/3oln(/xVA2) 

'Pi 



^ (3i In [ln(/xVA 



+ 



^2 In(MVA') Piln^f^yA^) 



{ln2(/xVA^)-ln[ln(/.VA^)]-l} + f 



(1.7) 



The dependence of the measurement of ^^(Mz) using Rz is illustrated in Fig. 2. The case 
of Rz nicely demonstrates the evolution of the scale dependence from NLO to NNLO. The 
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structure of the perturbative series indicates that a natural value for the scale /x is of order 
Mz, but its exact value is undetermined. 

In e+e~ annihilation it has become customary to take a central value of = 1 and 
to estimate an uncertainty in the prediction by varying Xfj_ between ^ and 2. (In contrast, 
in global analyses of parton distribution functions it is sometimes that is varied in that 
range, e.g. [9].) A variation of from ^ to 2 applied to i?z yields a theoretical uncertainty 
for as of ±0.0020, slightly lower than, but of the same general size as that obtained from 
the NNLO to NLO difference. This is despite the fact that the uncertainty from the 
variation is formally higher order in ag than that from the NNLO to NLO difference. 

In contrast to Ri-, event-shape distributions, which shall be studied in this paper, are 
known only to next-to-leading accuracy. In principle here too the last term of the series 
could be used to estimate the uncertainty. However, this is often seen as unduly conservative 
(it would give an uncertainty of order a^, whereas missing terms are of order a^) and 
the standard approach for NLO fixed-order calculations is to use just scale variations to 
estimate uncertainties. 

Also in contrast to the Rz case, over a significant part of the phase-space the event- 
shape distributions have divergent (or very poorly convergent) perturbative expansions, due 
to logarithmic enhancements of the perturbative coefficients. As a result, it is necessary to 
supplement fixed-order results with resummcd calculations [10-15], leading to a series with 
two expansion parameters — the coupling and the logarithm of the observable. Estimating 
uncertainties on the result is considerably more involved than for a traditional fixed-order 
expansion. For example, in the same way that the presence of a formally arbitrary scale 
Xfj, for the coupling is associated with an uncertainty, a new formally arbitrary scale xl 
appears in the definition of the logarithm and it too is associated with an uncertainty on 
the perturbative prediction. Further uncertainties arise from freedom in the way resummed 
and fixed-order predictions are combined ('matched'). 

A number of these possible sources of uncertainty have come to light only recently, 
partly in the context of DIS event-shape resummations [16]. There is therefore a need for 
a comprehensive up to date study of the whole range of such uncertainties for e+e" event 
shapes, as well as an understanding of how to combine them. That is the purpose of the 
present article. 

The definitions of the event-shape variables are given in Section 2, in Section 3 the 
theoretical predictions are summarised, in Section 4 the methods for estimating and com- 
bining theoretical uncertainties are presented, in Section 5 results obtained with the new 
method are discussed and the recommendation for an uncertainty strategy is summarised 
in Section 6. 

2. Definition of tiie event-siiape variables 

Event-shape variables are infrared and collinear safe observables that characterise global 
properties of hadronic events. They are experimentally measured with the momenta of 
charged and/or neutral particles. The variables are normalised so as to be dimensionless 
and independent of the production cross-section. One end of the spectrum is given by the 
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two-jet limit, the other by a kinematic hmit Umax, which in some cases corresponds to either 
spherical or symmetric event configTirations. The following variables are studied here. 

Thrust T [17]: the thrust axis unit vector Ht maximises the following quantity: 

T = n,^(S0^) . r=l-T, 

where the thrust sum extends over all particles i in the event. It is convenient to 
define r = 1 — T for the resummed expressions. 

Heavy Jet Mass p [18]: a plane perpendicular to ut divides the event into two hemi- 
spheres, Hi and H2, from which one obtains the corresponding normalised hemisphere 
(squared) invariant masses: 




The larger of the two hemisphere masses is called the heavy jet mass : 

P = M2 = -^max(M2,M|) , 

vis 

where Eyis is the total visible energy in the event. 

Wide Jet Broadening B\\r [19]: a measure of the broadening of particles in transverse 
momentum with respect to the thrust axis can be calculated for each hemisphere Hi 
using the relation: 

Hi- ^^-^ ,1 = 1,2 , 

where j runs over all of the particles in the event. The wide jet broadening is the 
larger of the two hemisphere broadenings: 

B]v = max(Si, B2) ■ 

Equivalently, the narrow jet broadening Bn is defined as B^ = min(Si, ^2). 

Total Jet Broadening Bt- the total jet broadening is the sum of the wide jet broaden- 
ing and the narrow jet broadening: 

Bt = By/ + B]\[ . 

C-parameter C [20] : the C-parameter is derived from the eigenvalues of the linearised 
momentum tensor : 

Ei \Pi\ V 1^*1 

The three eigenvalues Xj of this tensor define C with: 

C = 3 • (A1A2 + A2A3 + A3A1) . 
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Durham jet resolution parameter ys [21]: the Durham clustering algorithm for jet 
rates is taken as follows. For each pair of particles i and j in an event the metric yij 
is computed: 

2mm{Ef,E^){l- cos Oij) 
yij = (Durham algorithm) . 

VtS 

The pair of particles with the smallest value of y^j is replaced by a cluster. The four- 
momentum of the cluster is taken to be the sum of the four momenta of particles i 
and j, = Pi +Pj ('E' recombination scheme). The clustering procedure is repeated 
until all yij values exceed a given threshold ycut- The number of clusters remaining 
at this point is defined to be the number of jets. 

The ys = y2s jet resolution parameter is the threshold value of ycut below which 
an event is classified as having two jets and above which it has three jets. The 
distribution of ys falls steeply, with only a very narrow peak at small 2/3. In order to 
better examine the region of small ys, the logarithmic form — In 2/3 is usually analysed. 



3. Theoretical predictions 

In this section the ingredients of the theoretical calculations as well as methods for gauging 
the uncertainties are outlined. At the end of the section directions in which there is 
partial theoretical progress in improving the accuracies of the theoretical calculations are 
mentioned. 

3.1 Theoretical ingredients 

Fixed order calculation To second order in as, the distribution of a generic event-shape 
variable y {y=T, p, Bt, Bw, C or ys) is given by: 

J_^ = a,{^2^A{y) + {a,{,,')f[A{y)2Polnxl + B{y)] , (3.1) 
o'tot ay 

where cig = — , a;„ = , /x = renormalisation scale. (3-2) 

2tt Q 

The coefficient functions A and B are obtained from integration of the ERT matrix ele- 
ments, using for instance the integration program EVENT2 [22]. Consider the cumulative 
cross section: 

R{y,as) = / 1 -dx , (3.3) 

crtot Jo dx 

which may be cast into the second-order form 

(y, a,) = 1 + A{y)as + [A{y) 2/?o In ^ + Biy)] a^'' , (3.4) 

where A and B are integrated forms of A and B, and the explicit scale dependence of 
has been dropped for clarity. 
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Resummed calculations For small values of y, the fixed-order expansion, Eq. (3.1) 
fails to converge, because the fixed-order coefficients are enhanced by powers of Inl/y, 
A{y) ~ In^ y, B{y) ~ In^ y. To obtain reliable predictions in the region of y ^ 1 it is 
necessary to resum entire sets of logarithmic terms at all orders in a^. Certain event 
shapes and jet rates have the property that double logarithms exponentiate allowing one 
to write 

R{y, as) = (1 + Cias + Csa, + . . .) e^9i{a.L)+g2{asL)+a.9s{a.L)+... ^ ^^^^^^ ^ (3 5) 

where L = InyQ/y, with yo = 1 for y = t, p, y^, Bt, B\\r and yo = 6 for C-parameter. The 
function gi{asL) resums leading logarithms (LL), while g2{asL) resums next-to-leading 
logarithms (NLL), etc. Such a resummation scheme allows to make reliable predictions 
down to the region agL ~ 1.^ The gn{cXsL) functions have expansions 

00 

gn{asL) = J2 Gi,i+2-na/L'+^-'' • (3.6) 
1=1 

Expressions for the LL and NLL functions gi and 52 have been derived for a range of 
observables [10-15]. Additionally the Ci coefficients are known analytically, and the further 
subleading terms G21 and G2 are known numerically [10-15,23]. Thus the current state of 
knowledge for resummed results can be written 

RnlUv, as) = (1 + GiOs + C2as^) exp{Lgi{asL) + g2{asL) + G2iLas'^) . (3.7) 

The full set of coefficients to 0{a^) is given in Table 1. 





T 


P 


C 


Bw 


Bx 


-In y-i 




1.053 


1.053 


5.44 


1.826 


1.826 


-6.685 


C2 


34 


40 


76.5 


116.3 


92 


18.2 


Gn 


4 


4 


4 


8 


8 


4 


G12 


-2.67 


-2.67 


-2.67 


-5.33 


-5.33 


-4/3 


G21 


22 


36 


63.4 


73.8 


78.5 


-7.2 


G22 


-24.94 


-13.24 


-24.94 


-15.09 


-61.88 


0.868 


G23 


-10.22 


-10.22 


-10.22 


-27.26 


-27.26 


-3.407 



Table 1: Numerical values of the resummation coefficients, to 0{a^g) for n/ = 5. 

Matching fixed order to resummed calculations Pure fixed-order expansions are 
valid from moderate y to large y (a^. In^ y <C 1), while resummed calculations apply to 
small y (y <C 1). To obtain predictions over the whole kinematical range it is necessary to 
match the two calculations. This involves adding the two calculations and subtracting off 

^It is to be noted that two difTerent logarithmic classifications schemes are in use in the literature. The 
convention adopted here classifies the logarithms in Ini? and stems from [11]. In certain other contexts 

(e.g. [21]) it is logarithms in R itself that are classified, so that LL moans terms a^L^" in R, NLL gives 
terms a"I/^"~^, etc. Such a resummation scheme is valid in a more restricted range of L, asL^ < 1. 
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double counting. This alone is not sufficient because even after the subtraction of double 
counting, there remain terms from the 0{a^) contribution, in particular a piece G2i(x^L, 
which would cause the matched event-shape distributions to have a G2ia^/y divergence 
at small y. In contrast, the physical requirement is that the distribution should vanish at 
least as fast as a positive power of y. The matching procedure is therefore more involved 
than a simple subtraction of double-counting terms between the resummed and fixed-order 
contributions; the details are given below. 

Matching can be performed either for R or the logarithm of R, the resulting expressions 
are identical to 0{a1), but differ in the treatment of subleading terms. The prediction of 
the Log(R) matching scheme is given by [11]: 



lnR{y,as) = Lgi{asL) + g2{asL) - {GuL + Gi2L')as 



(3.8) 



- {G22L' + G2zL^)as'' + A{y)as + 



m - ^A\y) 



As the entire 0{a^) term, is exponentiated by this procedure, the problem of un- 

physical divergence from the G21 term is avoided. 

The expression for the R matching scheme reads [11] 



R{y, as) = (1 + Cias + C2Q!s^) exp [Lgi{asL) + g2{asL) + G2iLas'] 
- G21 La,2 - [Ci + GiiL + Gi2L2] a, 

C2 + Ci(GiiL + G12L') + ^(GiiL + ^12^2)2 + (G22L2 + ^23^3) 



(3.9) 



Us 



+ A{y)as + B{y)as'' . 

Here the G21 term is explicitly placed in the exponent, with the non-exponentiated fixed- 
order remainder vanishing as y is taken to zero. 

The Log(R) scheme is generally preferred over the R scheme, because the latter requires 
explicit knowledge of G21 and C2, which have to be evaluated numerically. The Log(R) 
expression is furthermore much simpler and it is believed to be theoretically more stable. 
This can be seen in the region of small y, where the remainder terms of the R scheme are 
found to be larger than the corresponding Log(R) terms. 



3.2 Sources of arbitrariness 

Modified matching The predictions obtained with the Log(R) and R matching schemes 
suffer from a limitation: unlike fixed-order predictions, they do not vanish at the multi-jet 
kinematic limit. The 0{al) expressions vanish at the phase space limit for four-parton 
production. The matching schemes can be modified to overcome this drawback. To do 
this, a kinematic constraint is imposed to guarantee that the prediction of the distribution 
vanishes at a given value ymax- This means for the modified Log(R) [11] 



lni?(yinax) = 



J_ da{y) 
(^tot dy 



y — Vmax 



dR 
dy 



= 



(3.10) 



y — 2/max 
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To fulfil this constraint L is replaced by 




(3.11) 



The power p is called the degree of modification and is usually chosen equal to unity. It 
determines how fast the distribution is damped at the kinematic limit. The nominal values 
of ymax are obtained for thrust, C-parameter and 7/3 on the basis of symmetry arguments. 
For the other variables matrix element calculations were carried out and compared to 
various parton shower simulations using PYTHIA [24], the results of which are given in 
Table 2. Ten million events were generated for the Monte Carlo simulation. Since the 
aim behind the modified matching formulae is to extend the distribution up to the true 
kinematic maximum of the observable, the maximum of all ymax determinations is taken 
as nominal j/max 

value. However, insofar as the prescription for this extension is quite 
arbitrary, for studies of the theoretical uncertainties it will also be instructive to examine, 
as an alternative for each observable, the lowest of the j/max values found in Table 2 (see 
Section 4). 



variable 


r 


P 


-In 2/3 




Bw 


C 


comment 


ymax 


0.5 




In 3 






1 


theoretical maximum 


ymax 


0.4225 


0.4175 


1.100 


0.4075 


0.325 


1 


ME 4-partons 


2/max 


0.428 


0.394 


1.102 


0.397 


0.307 


0.994 


PS partons 


I/max 


0.434 


0.383 


1.103 


0.396 


0.295 


0.995 


PS hadrons 


ymax 


0.5 


0.42 


In 3 


0.41 


0.33 


1 


nominal value 


^max 


0.42 


0.38 


In 3 


0.39 


0.29 


0.99 


alternate lower value 



Table 2: Values of j/max at which the distributions vanish. The 'ME 4-partons' entries have been 
determined from binned distributions and their accuracy is therefore limited by the resolution of 
the binning. *In the case of — lnj/3 it is the minimal value, corresponding to y^^max = 1/3. 



While for the modification of Log(R)-matching the replacement of L with L in Eq. (3.8) 
is sufficient to fulfil the constraints of Eq. (3.10) this is not true for the R-matching. To 
modify the R-matching in addition, the matching coefficients Gn and G12 become functions 
of y such that: 

^(ymax) = , G'ii(2/max) = , G'2l(?/max) = . (3.12) 

This is achieved with the following modification: 



L{y) = - In 
V 

Gniy) = Gil 
621(2/) = G21 



m 
LV y 

1 - 
1 - 



y 



2/max 

y 



2/0 

2/max 



+ 1 



(3.13) 



2/n 
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Finally the expression for the modified R matching scheme can be written as 



R{y, Us) = (1 + CiOs + C2as 



X exp 



Lgi{asL) + g2{asL) 



GiiagL + G2iLas^ 



G2iLas^ - 



C\ + GiiL + G12L 
1 



a. 



C2 + Ci{GuL + GuL^) + -{GuL + G^L^f + {G22L^ + ^23^^) 



+ A{y)as + B{y)as^ . 



(3.14) 



Renormalisation scale dependence For scale parameters different from unity, ev- 
ery second order terms acquires a scale dependence explicitly given by 



m 
G21 
G22 
C2 

92{oisL) 



B{y)=B{y) + 2^^A{y)\nxl, 
G2i = G2i + 2/3oGiilnx2 , 
G22 = G'22 + 2/3oGi2lnx2 , 
C2 = C2 + 2f3QCihixl , 

g2{oLsL) = g2{asL) + — {agLf g[{asL) In xl 



(3.15) 



The overlined terms of Eq. (3.15) replace for 7^ 1 the corresponding terms in equa- 
tions (3.8) and (3.9). Of course the coupling constant itself exhibits the scale dependence 
indicated in Eq. (1.7). 

Rescaling resummed logarithms In addition to the arbitrariness in the choice of 
there is also arbitrariness in the definition of the logarithms to be resummed; for example, 
whether powers of In ^ or of Og In ^ are resummed. This can be formalised [16] by 
the introduction of an xl parameter, analogous to Xn, such that where normally powers of 
as In — are resummed instead, powers of In are resummed. Such a rescaling alters 
the resummed formulae in the modified case according to : 



L 

92{0!sL) 



yo 



L = -ln 

P l\XL-y 

91 = 9i{oisL) , 



2/0 



•^L ' J/max 



+ 1 



92 



g2{asL) + InxL— {Lgi{asL)) 
dL ^ ' 



(3.16) 

(3.17) 

(3.18) 



Lgi{asL) I is referred to in some contexts as R'{asL). Rescaling 



where the quantity — ^ 
the argument of the logarithm also entails changes to the fixed-order coefficients both in 
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the modified and unmodified cases: 




Gl2 = Gi2 



(3.19) 



Gn = Gn + 2Gi2lnxL 



G23 = G23 



G22 = G22 + 3G23 In XL 

G21 = G21 + 2G22 In XL + 3G23 In^ XL 

di = Ci + Gn In XL + Gu In^ xl 

C2 = C2 + (CiGii + G21) InxL + (C1G12 + G22 + kGli) In' xl 



+(G23 + Gi2Gn) ln3 xl + ^Cfs In^ xl . 



Transformations of the expressions under x^ and xl variations are commutative and can 

therefore be carried out in any order. In the case of the modified R matching scheme, 



variation. 

3.3 Further estimates 

Direct estimates of higher-orders It is to be noted that for the thrust and heavy 

jet mass, investigations have been carried out of potential sources of higher order terms 
(NNLL, etc.) in the resummation, in particular those associated with the running of the 
coupling, using the dressed-gluon exponentiation model [23] . These estimated higher orders 
have a rather large effect on fits for a^, somewhat larger than the uncertainties which are 
deduced in Section 4. Given that these are strictly speaking only model calculations and 
that they exist for only a subset of the observables studied here (the thrust, heavy jet mass 
and C-parameter), they are not included here in the uncertainty estimates; their existence 
should however be kept in mind, together with the possibility that 'standard' methods (e.g. 
scale variations) for estimating the size of higher-order effects may be overly optimistic. 

Heavy-quark effects At Mz, events with primary b quarks represent about 20% of all 
events. The theoretical calculations assume however light quarks. It is important therefore 
to understand the impact of heavy quarks on the theoretical predictions. Fixed order 
predictions with heavy quarks have been in existence for a few years [25] and distributions 
with b quarks are known to differ by a few percent from light-quark distributions. Once 
this effect is multiplied by the fraction of b quark events it becomes of the order of a 
percent [26], which is small relative to the other perturbative uncertainties. 

Recently the first NLL resummed calculation for jet rates in heavy-quark events was 
completed [27] (though it is NLL for R as opposed to InR, i.e. a lesser accuracy than that 
used throughout this paper). Physically there are two effects at play. Firstly there is the 
'dead cone' (the suppression of collinear emissions with an angle smaller than rub/Q) which 
implies a modification of the double-logarithmic structure for below a critical value, y^c, 
of order ml/Q'^: 



factors of the type 1 — 




are to be applied to Gu and G21 (Eq. 3.13) after the xl 





(3.20) 
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Secondly higher order terms are modified because the number of active flavours (e.g. in 
/3q) decreases by one for the parts of the momentum integral with transverse momenta less 
than mb (corresponding also to < ysc)- The authors of [27] quote effects for jet rates of 
the order of a 3-4% for ycut = 0.004 with m;, = 5 GeV and Q = Mz- 

Currently no calculations exist for other observables. However the same physical ar- 
guments allow one to make the statement that the dead cone will lead to a modification of 
the double logarithms analogous to that of (3.20), below a critical value yc ~ m^/Q", with 
n = 1 for the broadenings and n = 2 for the thrust and C-parameter (the heavy jet mass 
is more complicated because of the direct b-quark mass contribution). In addition, for the 
thrust, jet-mass and C-parameter, in the range rrib/Q < y < ml/Q^, there is a mixture in 
the resummation of contributions with 4 and 5 active flavours. 

NNL calculations One source of future improvement in the theoretical accuracy is ex- 
pected to come from calculations of higher order contributions. Considerable progress has 
been achieved in calculating two- loop amplitudes for the process e+e~ qqg [28] . Subtrac- 
tion methods at NNLO to cancel infrared and collinear divergences between the two-loop, 
one- loop and tree- level contributions are currently being developed (see for example [29]). 
It will then be necessary to combine the various elements in the form of a fixed-order 
Monte Carlo program (analogous to EVENT2 [22]), which can be used to calculate the 
0{a^) contribution to the event-shape distributions. This is expected to reduce the scale 
dependence (see for example [30]), especially in the three-jet region. In the two-jet re- 
gion, the gain from the NNLO calculations could be more modest, because much of the 
contribution is already embodied in the NLL resummation. Improved accuracy over the 
full phase-space may therefore also require a NNLL resummed calculation. Progress on 
such calculations is also being made, though full results have so far been obtained only for 
observables that are somewhat simpler than event shapes, such as the Higgs pt distribution 
at hadron colliders [31]. 

4. Estimating theoretical uncertainties 

For the purpose of measurements it is necessary to adopt a nominal theoretical predic- 
tion, which is to be used to determine the central value of ag, as well as a set of variations 
for estimating the uncertainty on a^. 

For the matching scheme the use of modified Log(R) matching is advocated, since it 
has proven more stable than the modified R matching scheme, which will be used as the 
'alternative' theory for uncertainty estimates. 

The default value of ymax is the maximal possible value for any number of partons (as 
obtained from theoretical arguments and parton-shower simulations. Table 2), while the 
lower limit y'^^ax is used as an extreme alternative estimate. 

The options for the modification degree p are less clearly delimited — properties of 
certain fixed-order calculations [10,11,13] suggest that p should be > 1. The simplest case 
p = 1 is recommended for the nominal configuration. The effect of the value of p on the 
extracted value of ct^ is studied by fitting predictions with ag as free parameter and p ^ I 
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Figure 3: The change in extracted as as a function of the degree of modification -p for several 
observablcs. As reference theory the modified Log(R) matching scheme is taken with p = 1. 
Predictions with different p are fit to the reference with as as free parameter. 



to the reference theory with p = 1, using weights proportional to the differential cross 
section. The change in ag is depicted in Fig. 3. In the limit of large p the prediction turns 
into the disfavoured unmodified matching scheme. In general the largest difference in ag 
is observed around p = 2, which is suggested for systematic purposes. 

For the renormalisation scale, in a fixed order framework fi^ = is the simplest, 
hence natural choice for the invariant scale of the process. So in order to enable straight- 
forward comparison with existing measurements the conventional = 1 is recommended, 
together with the standard variation range, 1/2 < < 2. The renormalisation scale is 
kept the same in the fixed order and resummation parts of the calculation — though not 
strictly required this is also conventional, and varying separately in the two parts of the 
calculation would complicate the formalism somewhat. 

A novel type of systematic study is proposed for the proper resummation part of the 
theoretical prediction. The arbitrariness of the logarithmic terms to be resummed to all 



-14- 



orders is formalised through the variable rescaling factor x^. In analogy with the situation 
for renormalisation scale dependence, if a resummation to all orders of logarithmic accuracy 
(LL, NLL, NNLL, . . . ) was complete then the predictions would be independent of the 
choice oi XL- But at any truncated order (e.g. NLL) there is a residual xl dependence, 
which like the dependence, can be used to gauge the expected order of magnitude 
of missing higher order contributions. As in the case of the x^ scale dependence, the 
range of scales can not be derived from first principles. One of the most critical elements 
in the measurement of as and estimation of the associated uncertainties is the choice of 
a default value and range of variation for xl- All existing e+e" calculations implicitly 
use XL = i and this convention has always been assumed for measurements of a^. It is 
possible to argue for this as the most sensible default choice on the following grounds. In 
the resummation procedure, one evaluates integrals over the rapidities t] and transverse 
momenta of gluons. With the values of yo given after Eq. (3.5), one can show that 
the value of the observable in the presence of a single soft and collinear gluon emission 
is y ~ yo{kT/Q)^ • exp(— 677), where a and b are integers. Accordingly the logarithm that 
is resummed, L = hiyo/{xLy), can be rewritten alnQ/kt + brj — InxL- It then appears 
quite natural to choose the convention xl = 1, since L reduces to the combination of the 
physical logarithms, alnQ/kt + bf], without any extra constant piece. 

While there exists a simple motivation for the choice of a central value of xl = 1, 
implicitly embodied in current standard practice, the choice of the range of xl is far more 
subjective and has never been considered so far. An ad hoc prescription would be to 
vary xl in the same canonical range as X/^. The effect of a given value of x^ on the 
distribution must be studied thoroughly, however, the objective being to find an estimate 
of missing higher orders. A reasonable range for xl variations should not over-estimate 
the theoretical uncertainty, but complement the other investigations. Even if a strict range 
setting is impossible, a sensible proposal will be elaborated on the basis of various tests. 
This issue is studied in detail in the next subsection. 

Once all these different sources of uncertainty have been examined, it is necessary to 
combine them, bearing in mind that there is only partial complementarity between them. 
This is achieved with the uncertainty band method, which is discussed in Section 4.2. 

4.1 Setting a range for xl 

The determination of range of variation for xl is quite critical, because its effect is rather 
large. One approach is to try and find some theoretical motivation for a range, which leads 
to a number of possibilities: 

• The observation that many of the theoretical calculation involve an inverse Mellin 
transform that leads naturally to logarithms of e^^/y, suggesting a range | Inx^l < 7e, 
where Je ^ .57721566. This leads to 0.56 < xl < 1.78. 

• Certain values of xl lead some of the subleading fixed-order expansion coefficients 
being zero; for example Iuxl = | (twice this for — In 1/3) gives Gu = 0, suggesting 
|lnxL| < f, i.e. 0.47 < xl < 2.12 (double for -Inya). The appearance of a doubled 
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range for — log ^3 is natural also with the observation that for a single emission, the 
two-jet resolution parameter is the square of the jet broadening. 

• If instead G21 is wanted to be zero then a quadratic equation has to be solved for 
In XL, giving the results shown in Table 3. Taking the solution closer to zero as more 
natural suggests an average range of about | Inxil < 0.6, i.e. 0.54 < xl < 1.86. 



Observable 


(Inxi)- 


{lnxL) + 


T 


-1.98 


0.36 


C 


-2.39 


0.76 


Bj' 


-2.00 


0.48 


Bw 


-1.15 


0.78 


-In 2/3 






P 


-1.60 


0.73 



Table 3: Values of Ihxl which give G21 = 0. The results typically have an error of about ±0.05 
because G21 is known only from fits to fixed-order Monte Carlo results. In the case of — In 2/3 the 
solutions are complex. 

Another approach consists of a comparison of the 0(0^) calculation with the expansion 
to second order of the resummed NLL prediction. The latter is obtained by expanding 
Eq. (3.7), keeping only terms up to 0{a'^). The difference between these two expression 
is sensitive to asymptotic terms present in the exact O(a^) calculation but absent in the 
NLL expansion. This Ansatz conserves the information of the differential distribution, it 
can be expected that in general the uncertainties are not constant across the spectra. The 
difference A(y) between the NLL expansion and the exact 0{a1) calculation is determined 
for a central value of = 0.12. Then a theoretical variation is constructed by adding (resp. 
subtracting) A{y) to the reference prediction (i.e. using the modified Log(R) matching 
scheme with x/, = 1). Finally the reference theory is used with variable xl SlS free parameter 
to fit the variation. In practice, three different fit ranges, given in Table 4, are used to test 
the stability of the procedure. The first range (nominal) covers experimental fit ranges for 
as, the second range (2-jct) is restricted to the semi-inclusive region and the third range 
(3-jet) comprises multi-jet production. 



Observable 


fit range 1 (nominal) 


fit range 2 (2-jet) 


fit range 3 (3-jet) 


T 


0.70 - 0.97 


0.87 - 0.97 


0.70 - 0.83 


-lny3 


2.6 - 7.2 


5.0 - 7.2 


2.6-5.0 


P 


0.03 - 0.25 


0.03 - 0.12 


0.14-0.25 


Bw 


0.04 - 0.20 


0.04-0.12 


0.12 - 0.20 


C 


0.08 - 0.70 


0.08 - 0.30 


0.38 - 0.70 


Bx 


0.05 - 0.28 


0.05 - 0.14 


0.15 - 0.28 



Table 4: Fit ranges used for the determination of upper and lower bounds oi xl- 
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For the fitting procedure statistical weights scaling with the square root of the distri- 
bution value are applied. The resulting limits on xl are given in Table 5. It turns out that 
they strongly depend upon the range of the fit: a variation o{ xl leads to a change in shape 
which is very large below the peak region in the two-jet limit. The effect is minimal around 
the peak and then increases continuously towards Umax- The shape of A(y) is similar, 
but its slope is much steeper. This is reflected in the dependence of the results on the fit 
range. The resummation technique and the evaluation of A(y) arc generally considered to 
be applicable in the semi-inclusive region, substantiating results obtained with the two-jet 
fit range. As a cross-check the same procedure is applied to the case of x^. In average and 
for the two-jet fit range a span for from 0.4 to 2.9 is found, in reasonable agreement, 
although slightly over-estimating the canonical range from 0.5 to 2. 



Observable 


XL nominal range 


XL 2-jet range 


XL 3-jet range 


T 


0.55 - 1.89 


0.61 - 1.67 


0.49 - 8.1 


-In 2/3 


0.28 - 2.6 


0.29 - 1.71 


0.25 - 2.86 


P 


0.48-2.17 


0.59 - 1.69 


0.27-4.9 


Bw 


0.36-3.10 


0.38 - 2.70 


0.33 - 5.0 


C 


0.36 - 4.40 


0.38 - 3.78 


0.35 - 14.2 


Bx 


0.50 - 2.47 


0.57-1.99 


0.47-7.1 



Table 5: Ranges oixL obtained with fits to the difference between the fixed order 0{al) calculation 
and the second order expansion of the NLL prediction, for three different fit ranges. 

The comparison with the x^ variation can be investigated by determining xl values 
such that on average, for a given set of observablcs (and fit ranges), the impact of the xl 
variation on a fit of is the same as that of the conventional x^ variation. While at first 
sight this may seem to make the xl variation procedure redundant, it should really be 
considered as a method for setting a conventional variation range. Different observables 
may then have sensitivities to the common xl range that may be (and in practice often 
are) quite different. Furthermore, the y-shape of the xl and x^ variation will be different 
even if their average effect on a fit of ag is the same, because they probe different subsets 
of possible higher-order corrections. 

The procedure of this method is as follows: two predictions are calculated with the 
modified Log(R) scheme, = 0.12 and xl = 1, one with x^ = 0.5, the other with X/j, = 2. 
Each of them are then fitted with the same theory but x^ = 1 and xl being the free 
parameter. The results are given in Table 6. 

Having considered this variety of criteria for choosing the xl range, it is proposed to 
set a convention for the xl range of | < < | (equivalently |lnxL| < 0.405), which 
gives an average uncertainty on Ug which is similar in magnitude to that from x^ variation. 
This is a slightly narrower range than comes out from the purely theoretical arguments 
and from some of the other tests. For the purpose of a more conservative estimate of the 
uncertainty, or if one wishes to consider the uncertainty on the uncertainty, it is suggested 
to examine also a wider range |lnxL| < 0.6. Both ranges will be used for the numerical 
estimates in the next section. 



-17- 



W Uocl VdiUlt: 


'y T n ^~\Tvi 1 n o 1 von cfa 
X llUlllllid,! idllgc 




'T' I- -IT di" T* O n (Td 

X O-J C L 1 dllgc 


T 


'Uyj -L.I 1 


— 1 7S 
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-In 2/3 


0.48 - 2.23 


0.27 - 2.64 


0.59 - 2.17 




0.58 - 1.82 


0.56 - 1.86 


0.82 - 1.38 




0.67 - 1.51 


0.63 - 1.56 


0.79 - 1.35 


C 


0.61 - 1.71 


0.57-1.81 


0.65 - 1.54 




0.67 - 1.53 


0.65 - 1.54 


0.71 - 1.46 



Table 6: Bounds of xl obtained with a fit to the theoretical prediction with = 0.5 resp. 
Xi, = 2.0. 

4.2 Uncertainty band method 

The new method to assess the theoretical systematic uncertainty for the measurement of 
as, called hereafter uncertainty band method, is composed of two main building blocks: 

1. A nominal reference theory, the modified Log(R) matching scheme, used experimen- 
tally to determine the value of a^. 

2. A collection of theoretical uncertainties (variations of the theory) of the event-shape 
distributions, used to derive the perturbative uncertainty of as- 

The following variations of the theoretical predictions for the distributions are taken into 
account: 

• the renormalisation scale is varied between 0.5 and 2.0, 

• the logarithmic rescaling factor Xl varied in between 2/3 and 3/2 

(for — log(y3) an equivalent effect is obtained with squared endpoints, i.e. a variation 
from 4/9 to 9/4), 

• the modified Log(R) matching scheme is replaced by the modified R matching scheme, 

• the nominal value of the kinematic constraint j/max is replaced by the lower limit y^ax 
and 

• the first degree modification of the modified Log(R) matching scheme {p = 1) is 
replaced by a second degree modification {p = 2). 

The uncertainty band method gives a direct relation between the uncertainty of 
and the uncertainty of the theoretical prediction. Two pieces of information are required 
to calculate the systematic uncertainty: the measured value of the coupling constant, a^, 
and the fit range used for its extraction. With these elements in hand, the uncertainty can 
be computed without re-fitting the data. 

The method proceeds in three steps. First the reference perturbative prediction is 
calculated for the distribution in question using a given value for the strong coupling 
constant a^. Then all variants of the theory mentioned above arc calculated with the same 
value of a^. In each bin of the distribution, the largest upward and downward differences 
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with respect to the reference theory are collected. Hence the uncertainty is set by the 
extreme values of the theoretical variants. Variants which lead to similar but smaller 
effects are not double-counted. The largest differences define an uncertainty band around 
the reference theory. 

In the last step, the reference theory is used again, but with variable ag. In the spirit 
of this method, all valid theoretical predictions must lie within the uncertainty band for 
the fit range under consideration. Starting from the nominal a^, scans of ag are performed 
and the validity of resulting predictions is checked in each bin of the distribution. The 
largest and smallest allowed values of ag fulfilling the condition are used to finally set the 
perturbative systematic error. The method is illustrated for all variables in Figs. 4-9, with 
an input value of = 0.12. It might be argued that the condition for valid predictions 
lying strictly inside the uncertainty band is too tight, because the uncertainty is basically 
set by one single bin, which value is subject to statistical fluctuation of the numerically 
computed coefficient functions. An alternative operating mode of the uncertainty band 
method consists of a fit of the reference prediction to the uncertainty band envelope with 
as as free parameter. In this case weights have to be assigned to the bins inside the 
fit range. A convenient choice are statistical weights scaling with the square root of the 
distribution. 

As can be seen in Figs. 4 9, the total theoretical uncertainty, i.e. the uncertainty band, 
varies drastically across the distributions. The uncertainties are large, of the order of 30%, 
at both ends of the spectra and reach a minimum of a few percent in the peak region. 
In the range of experimental fits for the uncertainty of the distribution is between 5% 
and 10%. This turns into an uncertainty for between 3% and 6%. There is also a 
considerable spread among the observables, the uncertainties for the total jet broadening 
are in the central part twice as large as for — In 7/3. 

A closer look at the individual components of the theoretical error reveals two major 
contributions: the variation of the renormalisation scale and of the x^-scale. In the semi- 
inclusive region the variation generates the dominant uncertainty while in the three-jet 
region it is the variation oi x^. These two effects are clearly complementary. Going down 
in size of the uncertainty, the degree of modification as probed by the difference between 
p = 1 and p = 2 is important in the hard three-jet region, where this effect dominates. 
Both the kinematic constraint and the matching scheme uncertainty are rather small in 
the central part of the spectra and become somewhat larger at the multi-jet end. 

The uncertainty of ccg is derived from the range of values leading to 'valid' predictions 
inside the uncertainty band. The resulting uncertainty depends on the fit range, which 
usually encloses the central part of the spectrum. The extreme, but still valid reference 
predictions, shown as full lines in Figs. 4-9, are in general determined in the central re- 
gion, they touch the envelope of the uncertainty band at one point and remain inside the 
uncertainty band even outside the fit range. 
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Figure 4: Theoretical uncertainties for thrust. The grey area represents the total theoretical 
uncertainty for the distribution. The first curve indicates the reference theory with maximum resp. 
minimum as, which determines the theoretical uncertainty Aag for a given measurement a^. The 
other lines illustrate the different contributions to the theoretical uncertainty. 
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Figure 6: Theoretical uncertainties for heavy jet mass. The explanation is given in the caption of 
Fig. 4. 
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Figure 7: Theoretical uncertainties for wide jet broadening. The explanation is given in the 
caption of Fig. 4. 
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Figure 8: Theoretical uncertainties for C parameter. The explanation is given in the caption of 
Fig. 4. 
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Figure 9: Theoretical uncertainties for total jet broadening. The explanation is given in the 
caption of Fig. 4. 
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total uncertainty 


+0.0062 


+0.0028 


+0.0052 


+0.0066 


+0.0065 


+0.0077 


with 1 In XL 1 < 0.6 


-0.0059 


-0.0032 


-0.0052 


-0.0066 


-0.0061 


-0.0073 


total uncertainty 


+0.0062 


+0.0032 


+0.0052 


+0.0061 


+0.0066 


+0.0074 


willi 'III ' mclliod 


-0.0058 


-0.0041 


-0.0052 


-0.0057 


-0.0059 


-0.0064 



Table 7: Theoretical uncertainties for = 0.12 and fit ranges given in the first row. For each 
variable, the total theoretical uncertainty as well as uncertainties stemming from individual sources 
are given. In the last two rows concern results with more conservative theory assumptions. 



5. Results 

With the uncertainty band method as tool, all aspects of the theoretical systematic uncer- 
tainties for as can be studied in detail independent of a measured distribution. In Table 7 
results for the total uncertainty as well as for individual components are given for a repre- 
sentative fit range and = 0.12 as input. The total uncertainty range from 3% for — In 2/3 
to 6% for the total jet broadening, with an overall average of 5%. The uncertainties related 
to and XL are two-sided and in general asymmetric, the other uncertainties go in a single 
direction. The fact that the total uncertainty is larger than any of the individual ones is 
a consequence of the complementary collection of uncertainties of the distributions in the 
uncertainty band. A clear ranking in size of the different components of the theoretical Tin- 
certainty appears for all variables: variations of and xl are most important, followed by 
the matching scheme uncertainty, while the kinematic constraint and modification degree 
uncertainties are small and of similar size. Also given in Table 7 are uncertainty estimates 
obtained with more conservative assumptions, namely a larger range for xl (| Inx^l < 0.6) 
and a fit of the reference prediction to the uncertainty band envelope ('fit' method). 

Since the uncertainty is calculated with a fixed value of a^, the dependence of the 
result on the input value must be investigated. The input value, normally taken from an 
experimental measurement, depends on the observable and the centre-of-mass energy. For 
a theory defined up to O(a^), the uncertainty is formally expected to be at 0{a^g). The 
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matching of fixed order and resummed calculations, however, may alter the scaling with a^. 
The evolution of the uncertainty with as is shown in Fig. 10. In this case the symmetric 
uncertainty is analysed, i.e. the mean of the upward and downward uncertainty. A fit 
of the form a + ba^ gives a good description of the dependence of the uncertainty on Ug. 
The parameters a and b obtained with the fit are given in Table 8 for each observable. In 
general a is found to be substantially smaller than ba^, confirming that to a reasonable 
approximation the simple scaling discussed above is observed. 

The two main components of the theoretical uncer- 
tainty originate from the variations of and x^^. The 
size of the uncertainties depends crucially on the range 
of variation, nominally from i to 2 for Xu and from | to 



Observable 


a 


b 


T 


3.610-4 


2.99 


-ln?/3 


2.410^4 


1.51 


P 


2.710-4 


2.38 


Bw 


6.810-4 


2.71 


C 


5.810-4 


2.88 




9.710-4 


3.18 



Table 8: Parameterisation of the 

uncertainty evolution with ag ac- 
cording to a form a + ha'l . 



I for XL. The dependence on the range choice is stud- 
ied by re-scaling the variation range by a factor of xr. 
Hence, the re-scaled range has an upper limit of Xup * xr 
{xup = 2 resp. | for resp. xl) and a lower limit of 
Xdown/xR {xdown = 5 resp. | for resp. xl). The 
effect of a range variation is expected to be linear in the 
logarithm of the re-scaling factor xr and the results are 
given as function of InxR for both the upper bound (pos- 
itive error) and lower bound (negative error). The dependence of the x^ component of 
the uncertainty on the range of variation is shown in Fig. 11. The size of uncertainty 
increases rapidly with the width of the range, the positive error flattens at large ranges 
for — In 7/3, p and Bw- The shape of the dependence is similar for all variables. The case 
of the XL component is shown in Fig. 12. Here the rise at small variation ranges is even 
steeper than for x^. For the variable — Inys the evolution is significantly flatter than for 
all other variables, a similar effect would have been observed with a quadratic re-scaling of 
the variation range. 

Finally, the dependence of the uncertainty on the fit range experimentally used to 
extract as from the distribution is studied. For this purpose the uncertainty band method 
including all uncertainty sources is applied to each bin of the distribution, as shown in 
Fig. 13. In the central part of the distribution the uncertainty is reasonably constant 
and stable with respect to variations of one or two bins at each end. Approaching the 
peak of the distributions or the extreme three-jet region, however, induces rapidly growing 
uncertainties. This dependence should be kept in mind when selecting a fit range. 
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Figure 10: Evolution of the theoretical uncertainty with a^. The lines are the result of a param- 
eterisation scaling with a^. 
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Figure 11: The dependence of the positive and negative uncertainty related to the variation of the 
renormalisation scale x^. The width of the variation range is modified such that the lower endpoint 
is rescaled by 1/xr and the upper limit by xr. 
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Figure 12: The dependence of the positive and negative uncertainty related to the variation of the 
logarithmic re-scaling factor xl- The width of the variation range is modified such that the lower 
endpoint is rescaled by 1/xr and the upper limit by xr. 
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Figure 13: The dependence of the positive and negative uncertainty on the value of the observable. 
All components of the uncertainty are included. 
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The theory uncertainty studies described here are all carried out on purely perturbative 
predictions. One could instead envisage carrying out these studies after correction to 
hadron-level. Whether this actually makes a difference or not depends on the details 
of how the hadronisation is included, for example whether as a bin-by-bin multiplicative 
factor, a simple shift of the distribution, or a transfer matrix. However, since hadronisation 
corrections will have similar effects both on the various alternate theory curves and on the 
reference prediction, the net impact of hadronisation is expected to be rather small and 
hence it is simplest and most transparent to carry out the analysis at a purely perturbative 
level 

6. Conclusion 

In this paper, a new method is presented for the assessment of theoretical uncertainties 
in as- This method evaluates the systematic uncertainty of the parameter from the 
uncertainty of the prediction for the distributions from which it is extracted. After a 
comprehensive review of theoretical predictions for measurements of as from event-shape 
distributions in e"'"e~ annihilation, the uncertainties of such predictions are estimated with 
an uncertainty band method which incorporates several variations of the theory differing in 
subleading terms and includes a new test for re-scaling the resummed logarithmic variables. 
As the uncertainty band method can compute these perturbative uncertainties of as inde- 
pendently of a measured event-shape distribution, it is especially suited for an unbiased 
combination of several observables or experiments. 

The recommended method for computing the central result is to use the modified 
Log(R) matching scheme with Xf^ = XL = l,p=l and values for ymax given in Table 2. 
The assessment of the perturbative uncertainty consists of a variation of the renormalisation 
scale X/j, from 0.5 to 2, of the logarithmic re-scaling factor xl from | to |, of a replacement 
of the modified Log(R) matching scheme by the modified R matching scheme, of the degree 
of modification p = Ihy p = 2 and of the kinematic constraint ?/max by its alternative y[aax 
given in Table 2. The different uncertainties should be combined with the uncertainty 
band method. Proposals are made for more conservative uncertainty estimates; these are 
a larger range for the re-scaling factor | lnxx,| < 0.6 and an alternative operating mode of 
the uncertainty band method. 

To combine several measurements from different observables or experiments the per- 
turbative uncertainties of the individual measurements should be re-estimated with the 
above method using a common input value for ac(Afz) in the uncertainty band method. 
This input value should of course be consistent with the final result obtained iteratively. 

Following these instructions even existing measurements of as can be equipped with 
an up-to-date estimation of their perturbative uncertainty and then used in a consistent 
combination. 

Acknowledgements 

We wish to thank Stefano Catani for helpful suggestions and Giinther Dissertori, Einan 



-32- 



Gardi, Klaus Hamacher and Oliver Passion for useful conversation on this subject. 



References 

R.K. Ellis, D.A. Ross and A.E. Terrano, Nucl. Phys. B 178 (1981) 421 

Yu.L. Dokshitzer and B.R. Webber, Phys. Lett. B 352 (1995) 451; 
Phys. Lett. B 404 (1997) 321; 

Yu.L. Dokshitzer et al., J. High Energy Phys. 05 (1998) 003. 

J. Abreu et al., DELPHI collaboration, Phys. Lett. B 456 (1999) 322; 
M. Acciarri ct al., L3 collaboration Phys. Lett. B 489 (2000) 65; 
PA. Movilla et al., Eur. Phys. J. C 22 (2001) 1; 

J. Abdallah et al., DELPHI cohaboration, Eur. Phys. J. C 29 (2003) 285. 

PA. Movilla Fernandez, hep-ex/0209022, MPI-PH-2002-12 

LEPQCD working group, note in preparation 

K.G. Chetyrkin et al., Phys. Kept. 277 (1996) 189; 
E. Tournefier, LAL 98-67 

LEPEW working group and LEP collaborations, CERN-EP-2002-091, hep-ex/0212036 

K.G. Chetyrkin et al., Phys. Rev. Lett. 79 (1997) 2184. 

S. Chekanov et al, ZEUS Collaboration, Phys. Rev. D 67 (2003) 012007. 

S. Catani et al., Phys. Lett. B 263 (1991) 491; 
Phys. Lett. B 272 (1991) 368. 

S. Catani, L. Trentadue, G. Turnock and B. R. Webber, Nucl. Phys. B 407 (1993) 3 

G. Dissertori and M. Schmehing, Phys. Lett. B 361 (1995) 167. 

S. Catani et al., Phys. Lett. B 295 (1992) 269; 
Yu.L. Dokshitzer et al., J. High Energy Phys. 01 (1998) 11. 

S. Catani et al., Phys. Lett. B 427 (1998) 377. 

A. Banfi et al., J. High Energy Phys. 201 (2002) 18. 

M. Dasgupta and G. P. Salam, Eur. Phys. J. C 24 (2002) 213; M. Dasgupta and 
G. P. Salam, J. High Energy Phys. 0208 (2002) 032. 

S. Brandt et al., Phys. Lett. 12 (1964) 57; 
E. Farhi, Phys. Rev. Lett. 39 (1977) 1587. 

T. Chandramohan and L. ClaveUi, Nucl. Phys. B 184 (1981) 365; 
L. ClaveUi and D. Wyler, Phys. Lett. B 103 (1981) 383. 

P.E.L. Rakow and B.R. Webber, Nucl. Phys. B 191 (1981) 63. 

G. Parisi, Phys. Lett. B 74 (1978) 65; 
J.F. Donoghue et al., Phys. Rev. D 20 (1979) 2759. 

[21] S. Catani et al., Phys. Lett. B 269 (1991) 432; 

W.J. Stirling et al.. Proceedings of the Durham Workshop, J. Phys. G 17 (1991) 1567; 
N. Brown and W.J. Sthhng, Phys. Lett. B 252 (1990) 657; 
S. Bethke et al., Nucl. Phys. B 370 (1992) 310. 



-33- 



S. Catani and M. Seymour, Nud. Phys. B 485 (1997) 291. 

E. Gardi and J. Rathsman, Nud. Phys. B 609 (2001) 123; 
E. Gardi and J. Rathsman, Nud. Phys. B 638 (2002) 243; 

E. Gardi and L. Magnea, J. High Energy Phys. 0308 (2003) 30. 

T. Sjostrand et al., PYTHIA, Comput. Phys. Commun. 135 (2001) 238. 

G. Rodrigo, A. Santamaria and M. S. Bilenky, Phys. Rev. Lett. 79 (1997) 193; P. Nason and 
C. Oleari, Nud. Phys. B 521 (1998) 237; A. Brandenburg and P. Uwer, Nud. Phys. B 515 
(1998) 279. 

The ALEPH collaboration, Eur. Phys. J. C 18 (2000) 1. 

F. Krauss and G. Rodrigo, Phys. Lett. B 576 (2003) 135. 

L. W. Garland, T. Gchrmann, E. W. N. Glover, A. Koukoutsakis and E. Remiddi, Nud. 

Phys. B 642 (2002) 227; 

S. Moch, P. Uwer and S. Wcinzicrl, Acta Phys. Pol., Ser. B: 33 (2002) 2921. 

S. Weinzierl, J. High Energy Phys. 03 (2003) 062; 
J. High Energy Phys. 07 (2003) 052. 

E. W. N. Glover, Nud. Phys. Proc. Suppl. 116 (2003) 3. 

G. Bozzi, S. Catani, D. de Florian and M. Grazzini, Phys. Lett. B 564 (2003) 65. 



-34- 



